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Q.^', Abstract 

' The state with a spontaneous interlayer phase coherence in a graphene based bilayer quantum Hall system is studied. 
This state can be considered as a gas of superfluid electron-hole pairs with the components of the pair belonging to 
^ different layers. Superfluid flux of such pairs is equivalent to two electrical supercurrents in the layers. It is shown that 
1^ the state with the interlayer phase coherence emerges in the graphene system if a certain imbalance of the Landau 
level filling factors of the layers is created. We obtain the temperature of transition into the superfluid state, the 
■ maximum interlayer distance at which the phase coherence is possible, and the critical values of the supercurrent. The 
' advantages of use of graphene systems instead of GaAs heterostructures for the realization of the bilayer electron- hole 
superconductivity is discussed. 
PACS: 71.35.Ji, 73.21.-b, 73.63.-b 
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I INTRODUCTION 

^ , The experimental discovery of graphene [1| served as a stimulus for performing large-scale investigations of the 
electronic properties of conductors with a Dirac dispersion law for charge carriers. On the one hand such investigations 
, are certainly fundamental. On the other hand the uniqueness of the properties of graphene makes it an extremely 
interesting material from the standpoint of different practical applications. 

In the present article we shall examine the prospects for using graphene to obtain superconductivity of bound 

■ electron-hole pairs in a system of two two-dimensional conductors separated by a dielectric. The possibility of such 
^ ' superconductivity was suggested quite a long time ago 0, Q . The system considered was a sandwich consisting of a 

two-dimensional n-type semiconductor ~ dielectric - two-dimensional p-type semiconductor. Somewhat later it was 
determined 0, 0] that a similar phenomenon can be obtained in a semiconductor heterostructure with two quantum 
wells (quasi- two-dimensional electronic layers). If such a system is placed in the quantizing magnetic field oriented 
perpendicular to the layers (a quantum Hall effect regime is realized) and the total filling factor of the Landau levels 

■ is 1 (i^i + = 1, where Vi is the filling factor in a layer), pairing of electrons in one layer with holes (unfilled states in 
I the lowest Landau level) in the other layer occurs. These pairs are bosons; as temperature decreases, such a pair gas 
' can transition into a superfluid state. Because pair gas is two-dimensional the transition will occur by the Berezinskii- 
' Kosterlitz-Thouless mechanism. Since the electron and hole components of the pairs belong to different layers, the 

CNJ ' superfluid flow of the electrically neutral pairs is equivalent to two electric supercurrents flowing in opposite directions 
in neighboring layers. 

The analogy with a quantum ferromagnet is widely used to describe such a state (which is called a state with 
spontaneous interlayer phase coherence) [a, 0, S H [13] ■ This analogy makes it possible to study, within the framework 
of the same approach, the case of zero imbalance of the filling factors of the layers {vi — V2 = 1/2) as well as a 
situation with nonzero imbalance (yi — v,vi = \ — v,v^ 1/2), which arises, for example, because of the external 
}J] ■ electrostatic field of a gate. For a very large imbalance (y ^ 1) the system can be described as a rarefi ed g as of bosons 
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\ (magnetoexcitons) with a dipole-dipole interaction. Such an approach was developed in Refs. 11, 12, Ij, 14, 15). We 
note that experimental investigations of two-layer electron-hole superconductivity have been performed primarily on 
quantum Hall system [H, [B, M, [3, [53, [HI, [ll [2" 



The superconductivity of electron-hole pairs in a grap hene-dielectric-graphene system was studied in Refs. [24 



H, [H [23, [H, [H, [13]. In Ref. [3 [H, [23, [13] and [2g the possibility of such pairing was studied in the absence 
of a magnetic field perpendicular to the layers (i.e. not in the quantum Hall effect regime). In the absence of a 
magnetic field the condition for BCS pairing of electrons and holes is that their Fermi surfaces must coincide. For 
semiconductors with a quadratic dispersion law for the charge carriers the last condition requires that the electron and 
hole effective masses be the same. For a Dirac carrier spectrum (which obtains in graphene) the required condition 
on the Fermi surface is satisfied automatically. Neglecting screening the theory predicts a very high (hundreds of 
degrees) temperature for the transition into the superconducting state in such a system [23 . [29| . Taking account of 
screening the estimate for the transition temperature is much less optimistic and lies in the millikelvin range [27I [28| . 
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In a quantum Hall system the screening effects should not be so strong, which gives hope of attaining high tem- 
peratures for the transition into the superconducting state. In addition, the quantum Hall effect itself in graphene 
is observed at high temperatures (right up to room temperatures) [sij . The question of the superconductivity of a 
gas of magnetoexcitons in a two-dimensional graphene system was examined in Refs. 
Refs. 



28, 3i 



28| and [30|. The approach of 



can be use only for low magnetoexciton densities {v <^ 1). But the low-density case is not optimal from 
the standpoint of reaching high transition temperatures. More likely, the maximum critical temperature will obtain 
for half- filling of the Landau levels in each layer. The present article is devoted to analyzing this problem. 



STATE WITH INTERLAYER PHASE COHERENCE IN A BILAYER GRAPHENE SYSTEM 

The starting point for studying electron-hole pairing in a quantum Hall system is the Coulomb interaction Hamilto- 
nian written in the lowest (active) Landau level approximation. Consequently, we shall begin by obtaining the desired 
Hamiltonian for the graphene system. 

To describe the properties of the electronic subsystem of graphene in a magnetic field we shall follow the approach 
presented in detail in the review [s^l (see also Ref. [33|). Graphene possesses a honeycombed crystal structure, which 
can be represented as two simple triangular lattices A and B inserted into one another. The distance a between the 
nearest carbon atoms (which belong to different sublattices) is 1.42 A. In the tight-binding approximation and taking 
account of tunneling only between the nearest sites the Hamiltonian of the system has the form 

^ = E + ^■^)' (1) 

where a^, 6^ (oj, hi) are, respectively, operators creating (annihilating) electrons on the i-th sites of the sublattices 
A and B, i « 2.8 eV is the tunneling amplitude, and a is the spin index. 

Having written the Hamiltonian ([1]) in the momentum approximation it is easy to verify that the energy band is 
divided into two sub-bands which possess only two nonequivalent points of contiguity K = (27r/3a, 27r/3\/3a) and 
K' = (27r/3a, — 27r/3\/3a) (the x axis is directed along the line connecting any pair of nearest-neighbor sites). In 
undoped graphene (containing one free electron per site) the Fermi level passes through the point of contiguity of 
these sub-bands. Correspondingly, the sub-bands can be interpreted as electron and hole bands. In connection with 
such a structure of the spectrum the low-energy excitations in undoped and weakly doped graphene can be described 
by adding a pseudospin index a = ±1, corresponding to states whose quasimomenta lie close to K and K'. For 
low energies the Hamiltonian is diagonal with respect to the pseudospin indices. It is also assumed that there are 
no interactions which destroy the diagonality with respect to spin. The characteristic quantum numbers of such a 
Hamiltonian possess a definite spin and pseudospin. 

For such states the Schrodinger equation written in the coordinate representation has the form 

where vp = ?>ta/2h. For graphene the indicated parameter equals w 10® cm/ s. An electronic state is a spinor whose 
components correspond to the sublattices A and B. Switching in Eq. ^ to the momentum representation we find 
the spectrum in the form E — ±hvpk (k is the modulus of the wave vector, measured from the point K or K'). 
Evidently, vp is the velocity of the electrons on the Fermi surface. For fermions with such a spectrum this velocity 
does not depend on the carrier concentration. 

In a magnetic field perpendicular to graphene Schrodingers equation ([2]) assumes the form 



The solution of Eq. ^ gives the Landau energy levels in graphene: Eq — 0, E±n = ±{hvF/i)V2N, where i = 
y^hc/eB is the magnetic length and N = 1, 2, ... In Eq. ([3]) we neglected the Zeeman terms, since for reasonable 
values of the magnetic field the Zeeman splitting is much less than the spacing between the Landau levels. The 
eigenf unctions ^N.a,k corresponding to the zeroth and ±7V-th Landau levels have the following form: 
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*±w,+i(x,2/) = — — — — ===e { ' ], (5) 



where X — k£'^, and Hn{x) is a Hermite polynomial. The degeneracy of a Landau level is 4:S/{2tt£'^), where S ~ L^Ly 
is the area of the layer. The factor of 4 is due to the two- fold degeneracy with respect to both spin index and the 
index a. We introduce the filling factor — 27r^^n, where n is the concentration of filled states. A completely filled 
level corresponds to = 4. In undoped graphene the chemical potential corresponds to the zeroth Landau level, so 
that the negative values are completely filled, the positive values are empty, and the zeroth level is half-filled (its 
filling factor is 2). 

Let us consider a system of two graphene layers separated by a dielectric layer with thickness d and permittivity 
e. To avoid complications which are not of a fundamental nature we shall assume that the system is located inside a 
dielectric matrix with the same permittivity e. The Coulomb interaction Hamiltonian for such a system has the form 

= / d'rd'r%,A\r-r'mr)pAr')- (6) 



where Vi^i'{r) — jie^Jr^ + c?^(z — i'Y) is the Coulomb potential, 



is the electron density operator, i = 1, 2 is the index of the layer, and ^ and *I'i,Q.cr(r) are operators creating 
and annihilating at the point r an electron with prescribed spin and pseudospin. 

We shall represent the operators 'I'"'", ^ in terms of the operators ([4]) and (O creating and annihilating electrons 
in the quantum states, respectively. If the separation of the Landau levels is much greater than the Coulomb energy 
e^/(£^), then it is sufficient to retain in this expansion only a single active (partially filled) Landau level 

k k 

where A is the number of the active level, 'i'x a k two-component vector-matrices determined by Eqs. @ and ([5]), 
a k hermitian-conjugate matrices, and ^ ^ and ai^k,a,a are operators creating and annihilating electrons 

in a state with the corresponding quantum numbers in the level A. 

In this approximation the Hamiltonian ^ written in the Fourier representation has the form 

where Vi^i'{q) = {iire^ /eq) exp(— qdji — is the Fourier component of the Coulomb potential. For A = the Fourier 
component of the electron density operator is 



a,a k 

For A 7^ 



P»(q) = EE4fe+g„/2,a,.a^fc-9«/2,«.-exp (-iq^kf - 

rv rr h ^ ' 



(9) 



(q) = Z^Z^<fe+9„/2,a,aa-fc-<J„Aa,-exp( - —I ^ , (10) 



a, a k 



where L\{x) is a Laguerre polynomial. 

If the system is composed of two undoped graphene layers, then in each layer the half-filled zeroth Landau level 
will be active. If an electrostatic field is applied in a direction perpendicular to the layers, then an imbalance of the 
filling factor can arise, i.e. the filling factors of the layers will become equal to z^i = 2-\-v and 1^2 = 2 — z> (0 < i/ < 2). 
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In stronger fields the active levels will become A = +1 in layer 1 and A = — 1 in the layer 2, then +2 and —2, and so 
on. 

We shall examine first the case where the zeroth levels are the active levels. The Coulomb interaction Hamiltonian 
diff'ers from the Hamiltonian of the usual bilayer quantum Hall system with total filling factor v ~ 1 only by the fact 
that it contains an interaction of not one but four components, corresponding to the quantum numbers (a = ±1, a =t 
, i). We shall enumerate these sets of quantum numbers by the index P — 1,2, 3, 4. 

By analogy with the single-component case 0] we write the trial multi-particle wave function of the system in 
the form 



+ 



k 13 



|0). 



The function (|lip describes a state where interlayer phase coherence arises for each component, 
(see, for example, Refs. 



(11) 



It is easily shown 



3J,|35[) that the state pTjl can be represented in the form of a BCS wave function describing 



the pairing of electrons of one layer with holes of the other layer with the paired electrons and holes corresponding 
to the same state 13. The parameter Op is related with the magnitude of the imbalance for this component by the 
relation zJ^ = cos0^/2; the filling factors for this component are i^i(2)/3 = 1/2 ± i/^, and the total imbalance is 

The energy of the system in the state (fill) is given by the expression E = {'^\Hc\^) + -Eg + i?bg, where is the 
interaction energy of the electrons with the electrostatic field produced by the external gate, and i?bg takes account 
of the Coulomb interaction with the positive core. A direct calculation (see, for example, Ref. [s^) gives 



E = 




- Jn 1 



E 

P 



-Ji 



(12) 



where erfc(a;) is the complementary error function and V is the potential difference produced between the layers by 
the gate field. 

The first term in Eq. (|12p is identical to the energy of a plane capacitor and represents the contribution of the 
direct Coulomb interaction to the energy. The second and third terms contribute to the exchange interaction inside 
and be tween the layers, respectively. We note that Jq > Ji, i.e. the in- layer exchange energy constant is greater than 
the interlayer exchange energy constant. 

To find the parameters 9/3 appearing in the trial function (|lip it is necessary to find Dp corresponding to the 
minimum of the energy (|12p taking account of the constraint \i>fj\ < 1/2 (which follows from the definition oi i>p). 

In a nongraphene bilayer quantum Hall system with I'l + = 1 there is only one component, and the energy 
minimum in the absence of a gate field corresponds to = 0, i.e. 1/1 = 1/2 — 1/2. The situation is more complicated 
in a multi-component system. For Jq > Ji the part quadratic in i>p in Eq. (|12p is not positive-definite. Moreover, 
it is not positive-definite with one and two fixed. Therefore the minimum of the energy (|12p is attained at the 
boundary of the range of i/p, and in addition at least three of the four i>p must assume limiting values, equal ±1/2. 
For V = the minimum is reached for 

E^/' = 0' 

but when two i/p- = -1-1/2 and the remaining other two vp., = —1/2. Physically, this result is understandable. Since 
Jo > Ji, a gain in the energy exchange is attained for maximum imbalance of a given component. In a one-component 
system a direct exchange interaction impedes this. In a multi-component system the latter effect is absent, since 
the magnitudes of the imbalance of the components can be of opposite sign. The values i>p = ±1/2 correspond to 
6p = 0,TT which means absence of interlayer phase coherence (i.e. absence of a superconducting state). We note that 
this situation is analogous to the one arising in a nongraphene bilayer quantum Hall system with total filling factor 
i/i + 1/2 = 2 [33|. 
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We shall show that for y ^ it becomes possible for a state with interlayer phase coherence to appear. For 
definiteness, we shall assume that V > 0. We shall also assume that the parameters determined by Eqs. (fT3)) satisfy 
the inequality W — Jq + Ji > 0. Analysis of Eq. shows that for eV < Jq — Ji an energy minimum corresponds 
to the same state as for V = 0. If the potential difference satisfies the condition 

Jo-Ji<eV<2W-Jo + Ji, (14) 

then the minimum energy is attained in the state 

The state ((T5|) is degenerate with respect to transposition of the indices (i. Most likely, such degeneracy is removed 
as a result of weaker interactions which are not taken into account in Eq. (|12p . This is not essential for our analysis. 
In the parameter range 

2W - Jo + Ji < eV <2W + ,h - Ji (16) 

the minimum corresponds to the state 

11 

Vl=V2=VZ^ -, = --. (17) 

If V lies in the range 

2W + Jo - Ji < eV < AW - Jo + Ji, (18) 

then the minimum is reached when 

I ^ eV-iW 
"^^"^ = "^=2' "^=2(W^-Jo + JO- ^ ^ 

As V increases further, the zeroth Landau level will become inactive: it will be completely filled in layer 1 and empty 
in layer 2. 

In the next section we shall see that the superfluid density and, correspondingly, the transition temperature are 

proportional to i^^(l — vp) — y^l/4 — v^^ i.e. different from zero for ^ ±1/2. As follows from the expressions 

presented above, pairing will arise when the external potential difference lies in the range (|14p and (|18p . Then the 
maximum superfluid density (corresponding to Vfj — 0) will correspond to eV = W and eV = iW . To estimate the 
required field strength in the gate we note that the condition eV = W corresponds to field strength E m 2 ■ lO^e^^i? 
V/cm, where B is taken in T. 

We also note that for eV = W the filling factors of the layers are i/i = 5/2 and 1^2 = 3/2, and for eV — 3W they 
assume the values vi = 7/2 and 1^2 = 1/2. In both cases only one component has a nonzero imbalance between the 
layers, and it is this (active) component that is responsible for the pairing of electrons and holes. 

As V increases further, the electrons will transition from the level A = —1 of layer 2 to the level A — +1 of layer 1, 
and pairing of electrons in the A = +1 level with holes in the A = — 1 level will become possible. According to Eq. 
(fTOl) . the expressions for the Fourier components of the density operator for the levels A = ±1 have the form 

a,(j k \ / \ / 

The state with pairing of electrons in level A = +1 with holes in level A = — 1 can be described by a wave function 
of the form (fTT|) in which the operators creating electrons in the layer 1 belong to the Landau level A = +1 and the 
operators creating electrons in the layer 2 belong to the level A = — 1. 

The Coulomb interaction energy in such a state, taking account of the imbalance of the components of the level 
A = 0, assumes the form 

E^A.iwi2 + ^(l + DS\ -2Jo-J2fl + ^z>^) 
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where 



-Ja 1 




2eV - {eV - 



(21) 



16 V 2 rf' 



J3 = J2 



11 - 2d^ 
IT" 



■ exp 




(22) 



are the in-layer and interlayer exchange interaction constants for the levels |A| — 1 {d — d/£). As in the case of an 
active zeroth level, the in-layer exchange constant J2 is greater than the interlayer exchange constant J3. In Eq. (PT|) 
O = 2^/2hvF/i is the splitting between the levels A = -1-1 and A = — 1, 7 denote the same quantum numbers as f3 but 
for the levels |A| = 1, and i/^ determine the imbalance of the components of 7. 
We find from Eq. ^ that for 



4W + n + J2-J3<eV<6W + n-J2 + J3 



the minimum energy corresponds to the state 

eV-n-5W 



2(M^- J2 + J3)' 



V2 = ^3 



1^4 



1 

2' 



(23) 



(24) 



i.e. the inequalities (j23p determine the next range of V where a state with interlayer phase coherence can arise but 
this time between the electrons belonging to the levels A = +1 and A = —1. Specifically, for eV — + 5W the 
component in which such coherence will arise will have zero imbalance, i.e. the maximum superfiuid density will be 
reached. 



CRITICAL PARAMETERS FOR ELECTRON-HOLE SUPERCONDUCTIVITY IN GRAPHENE 

As shown in the preceding section, a state with interlayer phase coherence can arise in a bilayer graphene system 
in a quantizing magnetic field and an electric field perpendicular to the layers. We established that even though 
four components are present in a Landau level in graphene, only one component can be active. Depending on the 
magnitude of the gate field interlayer phase coherence can arise between the electrons in the zeroth Landau level or 
between the electrons in the A — ±1 levels. In the first case the critical parameters of electron-hole superconductivity 
in graphene will be described by the same equations as in bilayer quantum Hall systems in GaAs heterostructures. 
In the second case the equations will be somewhat different. In the present section we shall cornpare two such cases. 
We shall use for the analysis the approach of Ref. [l3| , which was further elaborated in Ref . [s^ , as well as in the 
present work. 

We shall examine a state where V gives interlayer phase coherence of one of the components, and in which the 
parameters 9 and (p for this component depend on k: 

I*) = n ('^"^ Y«tfc + Bin y 4..) |0). (25) 

Here and below the creation and annihilation operators as well as the functions 9^ and (/s^ refer to the active component. 
The quantum index k determines the coordinate x of the center of the electron orbit: X = kl'^ . Consequently, the 
function describes a state in which the order parameter {'^\a^y.a2x\^) = (1/2) sin Oxe^'^^ for electron-hole pairing 
varies along the x axis. 

The energy of the system in the state (|25|) is 

E = Eo-eV^cos9x + ^ ^ ^[H{X - X') - F<i{X - X')]coBex cos9x' 

X y x.x' 

-FD(X-X')sin0xsin0x'cos(^jc -^x')}: (26) 

where V = V —Vq\s the potential difference corresponding to zero imbalance of this component and Eq is the Coulomb 
energy of the inactive components. To abbreviate the notations appearing in Eq. (|25p we shall express the functions 
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oiX-X' in terms of the relation A{X) = J dqA{q)e^'^^ , where A{X) = H{X), Fs{X) and Fd(X) and the exphcit 
form the Fourier component of these quantities is as follows: 

p2 2,2 1 _ p-d\q\ 

H(.) = ^e^-^^/.M), (27) 



dk t 



■Mkq£)h{k), 



(28) 



(g) - / dke--Jo (kqi) e-'^fxik) 



(29) 



In Eqs. ((27 |) -([29 | Jo('i') is the Bessel function of order zero, and 



hik) 



for A = 
for lAI = 1 



(30) 



(the active levels are A = OorA = ±l). 

In the uniform sate with nonzero flux of electron-hole pairs in the x direction the phase of the order parameter is 
linear in X {(px = QX) and 9x = Oq is independent of X. The energy of this state is 



E, 



mf ■ 



En 



S 
4^ 



[n{0) - cos^ 00 ~ ^d(Q) sin2 00 



(31) 



Since the imbalance vp = cos 6*0/2 and 27i:(0) = W, 2J^s(0) = Jo (J2), and 2J^d(0) Ji (J3) (for A = (±1) 
respectively), it is evident that there is complete correspondence between the expression ([31]) at Q = and the 
expressions and (PTjl . 

For small Q the energy (PT]) can be represented as 



where 



= 5(const + -paoQ^) 



,2/1 Ec fd 



Pso = sm c/o-— p - 
Idtt 



= e^/ei is the Coulomb energy, and the function 



(32) 



(33) 



p{x) 



yfe^erfc(^) (l + a:^)-^:, 



y|e^2 erfc(^^^ 



7+13x^+7a:*+a:'' 
16 



16 



for A = 
for |A| = 1 



(34) 



gives the dependence of this quantity on the distance between the layers. 

The quantity p^o is called the superfluid stiffness. More precisely, the expression (1331) gives the value of this quantity 
at T = in the mean- field approximation. A plot of the function pso{d/£) is displayed in Fig. [TJ The quantity ttpsq/2 
gives an estimate of the temperature of the transition into the superfluid state. This estimate neglects the temperature 
correction to ps and is valid in a bounded interval of d. We shall discuss this question in greater detail below. 

To find the critical current and the transition temperature it is necessary to obtain the spectrum of the collective 
modes. We shall use the approach of Refs. [l3| and [1^, which is based on the quantization of the energy of small 
fluctuations of 9x and (px ■ First we shall require the spectrum for wave vectors q directed along the gradient of the 
phase if (which was obtained in Ref. [11]). The general case is analyzed in the Appendix. According to Ref. jsl], the 
expression for the spectrum can be written in a form that is formally identical to the expression for the Bogolyubov 
spectrum of quasiparticles in a moving condensate: 



E{q) = ^/e{e + 2J) + hqv. 



(35) 
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FIG. 1: Dependence of the superfluid stiffness (in the units EcSin^ 60) on the distance between the layers (d = d/£). Solid line 
- A = (pairing on the zeroth Landau level), dashed line - |A| = 1 (pairing on the ±1 levels). 



The following notations have been introduced in Eq. ((35|) 



7 = sm Wo 



■H{q) - Ts{q) 



:FD{q + Q)+TD{q-Q) 



(36) 
(37) 



J'pjq + Q) - Tpjq - Q) 
hq 

The quantity p6p is the kinetic energy of pairs. In the limits q — > 
e = fi^q'^ /2M, where M is the magnetic mass of a pair 



cos Oq. 



(38) 



M 



2h^e 1 
p{d/iy 



and (5^0 the expression ([M]) reduces to 



(39) 



and p{x) is determined by Eq. (|34|) . Since the mass M is inversely proportional to ps, according to Fig. 1 it increases 
with d, and the magnetic mass for the active levels A = ±1 is greater than that for the active level A = 0. 
The expression (l37l) in the limits q ^ and Q — > goes to a constant 



70 



where 



7(2;) 



e 2 erfc-^ 



16 



/ 11 _ ll-2x^+^" 4 f 
\^16 16 ^ ^^^^ 



V2 



for A = 
for |A| 1 



(40) 



(41) 



The quantity (|40|) determines the pair interaction energy per pair. The function jo{d/£) is presented in Fig. [2l 
Evidently, the interaction energy for A = ±1 is greater than that for A = 0. Moreover, for small d the energy 70 cx; d 
for A = ±1, while for A = the function jo{d) is quadratic. 

In the long- wavelength limit the spectrum is linear: E = hsq ((5 = 0), where the sound velocity equals 



s = Vto/A/ = — sin 6I0S f - ) , 



(42) 



s{x) = -Vl{x)pix)- 

A plot of the function s{d/£) is presented in Fig. [3l According to Fig. [Sj the maximum sound velocity (which is 
reached at d/£ « 1 and zero imbalance of the active component) does not depend on the magnetic field and for e w 3.9 
(SiOa) it is « 10^ cm/ s. 
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FIG. 2: Average interaction energy between pairs (in the units Ec sin^ 60) versus the distance between the layers. Sohd hne 
A = 0, dashed line - |A| — 1. 




FIG. 3: Velocity of the acoustic mode (in the units /{eh) sin6'o) versus the distance between the layers. Solid line - A = 0; 
dashed line - |A| — 1. 



The expression ([55)1 in the limits q 0, Q ^ assumes the form v — {hQ/AI) cosOq. Since Q ~ dtp/dX, the 
quantity v is the product of the superfluid velocity (vs — hVip/M) and an additional factor that depends on the 
imbalance. Specifically, for zero imbalance v — 0. As discussed in Ref. [38], this is a consequence of the electron- hole 
symmetry. In the present work we want to elucidate this feature on the basis of several other arguments. 

The superfluidity of electron-hole pairs in a bilayer quantum Hall system can be viewed as an analog of the so-called 
counterflow superfluidity [39*1 . Counterflow superfluidity arises in a system in which type I and II bosons on a lattice 
are present, the total filling factor of the lattice sites is 1, and strong same-site repulsion forbids two bosons from 
occupying the same site. An elementary act of boson motion in such a system is a type-I boson hopping onto a site 
with a type-II boson and the type-II boson simultaneously hopping onto the site occupied by the type-I boson. The 
fluxes of type-I and -II bosons in such a system are equal in modulus and are oppositely directed. The ratio of the 
velocity of the components I and II depends on the concentration ratios of the components. For equal concentrations 
the moduli of the velocities are likewise the same. 

In the bilayer quantum Hall system considered here the number of electrons of the top layer is equal to the number 
of holes in the bottom layer, and likewise the number of holes in the top layer is equal to the number of electrons 
in the bottom layer. The pairing of electrons and holes from neighboring layers can be regarded as the formation of 
two kinds of pairs differing by the direction of the dipole moment (upwards or downwards). In such an interpretation 
pair motion is an exchange process between different kinds of pairs. In other words, the flux of pairs of one kind is 
accompanied by a counterflux of pairs of the other kind. In the absence of imbalance the concentrations of different 
kinds of pairs is the same. 

Counterflow superfluidity is a particular case of a two-component superfluid 4^ 41 , [43 | . In this particular case 
additional conditions are imposed-equality of the counterfluxes and conservation of the locally total density of the 
components (for this reason there is only one and two oscillation modes) . The spectrum of the collective modes in the 
two-component system possesses the feature that it does not contain terms which are linear in the gradient 

of the phase, if the velocities of the components are equal in magnitude and opposite in direction. At the same time 
this spectrum contains a quadratic dependence on the gradient of the phase. The spectrum (|35p demonstrates similar 
properties, and vanishing of v with zero imbalance is a consequence of the counterflow character of the superfluidity 
in a bilayer quantum Hall system. 

For finite values of q the spectrum deviates considerably from the Bogolyubov spectrum. Specifically, for sufficiently 
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FIG. 4: Energy of the collective mode (in the units Ec) versus the wave vector q (in the units £ ^) for 0/3 — and d — 0.9dc 
(!) and d = 0.99dc (2). Solid lines - A = 0, dashed lines - |A| = 1. 




FIG. 5: Critical distance between layers versus the imbalance of the filling factor of the active component. Solid line - A = 0, 
dashed line - |A| = 1. 



large d a roton-like minimum appears in the spectrum; the depth of this minimum increases with increasing d (Fig. 
|4]). For d > dc the spectrum becomes imaginary for finite q and an instabiUty relative to the formation of a charge 
density wave arises in the system. The dependence of the critical distance dc between the layers on the imbalance 
of the filling factors of the active component is presented in Fig. [5l whence it is evident that the component with 
pairing of the electrons and holes on the Landau level A = ±1 is unstable for smaller d than the same state arising in 
a situation where the active level is A = 0. 

The critical values of the supercurrents can likewise be determined by analyzing the spectrum ((35l) . As already 
discussed in Ref. [sl], the limitation on the currents in the present case is not associated with the magnetic field 
which the currents in the layers generate and is oriented parallel to the layers. The maximum value of the current 
corresponds to the maximum value of Q for which the spectrum psp is positive (Landaus superfluidity criterion) 
and real (stability criterion). The dependence of the densities of the supercurrents in the layers on Q is given by the 
expression [ssj 

el . 2^ dT-D{Q) 

The computational results for the critical current are presented in Fig. [6l The maximum critical current is proportional 
to the quantizing magnetic fleld. For B = \ T and e = 3.9 the maximum critical current density equals approximately 
« 1 A/m. 

In closing this section we shall estimate the dependence of the temperature of the transition into the superfluid state 
on d. Since the system considered here is two-dimensional, the transition into the superfluid state is a Berezinskii- 
Kosterlitz-Thouless transition. The critical transition temperature is determined by the equation Tc = ■Kps{Tc) /2, 
where ps(T)/2 is the superfluid stiffness at finite temperature. To find PsiT) we note that the superfluid flux density 
is is related with the gradient of the phase Q by the relation js = Ps{T)Q/h (for small Q). On the other hand the 
flux density can be found from the relation js = {l/Sh)dF/dQ, where 

F = i?,„f + i?,p + r^ln('l-cxp('-:^')') - (44) 
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FIG. 6: Critical current density (in the units /hel^) versus the distance between the layers for different magnitudes of the 
imbalance. Solid line - i> = Q\ dashed line - i> = 0.17; dot-dash line - v = 0.25; dashed line - v = 0.4 for A = (a) and |A| = 1 
(b). 



is the free energy. Here _Bzp is the energy of zero-point vibrations. The collective and single-particle excitations must 
be taken into account in the energy of the zero-point vibrations. As shown in Ref. ^8J, for d not too close to the 
critical value the contribution of zero-point vibrations to the renormalization of the superfluid stiffness is negligible. 
Consequently, we shall neglect i?zp when calculating the superfluid stiffness. To take account of the temperature 
correction to ps it is sufflcient to include only the collective modes in the entropy term in Eq. ()44p (the single-particle 
excitation spectrum has a gap of the order of E^. As a result we have 

P.(T) = P.o4li^„^E^^-(^(^))^ (45) 

where Nb{E) — (exp(i?/r) — 1)^^ is the Bose distribution function. For a single-component superfluid system with 
the dispersion law i?(q) = EQ{q) + tiq^Vs (where Vg — hQ/M is the superfluid velocity, and EQ{q) is the spectrum as 
Vs = 0) the expression ([IS]) reduces to the standard form Ps{T) = Pso — Pn, where 



Af5^ M V dE 
q 



E=Eo(q) 



is the normal density. In our case the expression for the normal density will be different, since the superfluidity is of 
a counterflow character. 

Here we shall focus on the calculation of the transition temperature with zero imbalance. In this case the spectrum 
for small Q can be written in the form 

E[ci) ^ Eo{q) + \a{c^)Q^ + . . . 

which gives 

p,(T)=p,o + ^Ea(q)iVB(£;o(g)). (47) 
q 

To calculate the coefficient a(q) in the expansion we shall require the spectrum of excitations for an arbitrary 
direction of the wave vector. A method for calculating the desired spectrum, generalizing the approach used in Ref. 
[sst . is described in the Appendix. According to the results obtained in the Appendix, for do = 7r/2 the spectrum 
assumes the form 



E{ci) = ^^c,,Q{e^.Q + 2^c,,Q), (48) 

where 

eq,Q = 2Td{Q) - TT,{\q,+ Qx\) - .Fodq - Qx\), (49) 



7q,Q 



H(q, Q) ~ Ts{q) + ^D(lq + Q^I) + ^D(lq-Q^I) 



(50) 
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FIG. 7: Critical temperature versus the distance between the layers. The solid curves 1 and 2 correspond to A = and |A| = 1, 
respectively. The dashed curves show the critical temperature in the mean-field approximation (neglecting the temperature 
renormalization of the superfluid stiff'ness). 



the functions Ts{q) and J-uil) are determined by Eq. (|28p and (|29p . x is a unit vector in the direction x, and 

n,t n, ^2 _4il-e-;^^cosM^ 

^(q,Q) = ^e 2 _ a f^(q£)_ (51) 

Using Eqs. (I48ll5ip to calculate the function a(q) we arrive are the following expression for Ps{T): 

Ps{T) = Pso + [ (2(7, + e,)T';,[0) (7, + I) (^;(^) + + e,V[q^ , (52) 

where 

I?(9) = ^e-^g^2e-'^VA(g^), 

= eq and 7^ = 7q.o- The dependence of the critical temperature on the distance between the layers, as calculated 
from Eq. is presented in Fig. [T] 

According to the results obtained, for intermediate values of d the expression Tc = irpso/l is a good approximation 
for the transition temperature; it is important to take account of the temperature correction for small d/t and for 
d close to the critical values. The maximum critical temperature is attained for intermediate values of d. The 
dependence Tc{d/£) is qualitatively similar to the dependence of the critical temperature on d/£. In a state where 
electrons and holes from the zeroth Landau level form pairs the critical temperature is higher than in a state where 
electrons in the level A = +1 form pairs with holes in the level A = — 1. 

We note that the results of the present section, which were obtained for the zeroth level, also describe the situation 
in bilayer quantum Hall systems in semiconductor heterostructures based on GaAs (and, correspondingly, partially 
repeat the results of Refs. [1, 0, B 0, HI ) . 



DISCUSSION 



The essential result following from the preceding section is that the maximum parameters of electron-hole super- 
conductivity in a bilayer graphene system in a quantizing magnetic field can be attained if interlayer phase coherence 
arises between electrons belonging to the zeroth Landau level. In this connection there arises the question of whether 
or not there are any advantages to using graphene instead of GaAs heterostructures. 

To answer this question we recall the general conditions under which a superfluid state of bound electron-hole 
pairs can be expected to appear in a bilayer quantum Hall system. Three characteristic energies can be identified 
in the present problem: the splitting loq between the Landau levels, the Coulomb energy Ec, and the interlayer 
tunneling amplitude t. The analysis performed in the present work is based on the assumption that the following 
strong inequalities hold: t <^ Ec and Ec loq. 

The first inequality make is possible to neglect tunneling when analyzing collective modes. In addition, even weak 
tunneling results in the formation of a system of vortices which are similar to Josephson vortices. If a bilayer system 
is used to transfer current from source to load (this is the configuration realized in Refs. 2^, 21, i^]), the vortices 
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FIG. 8: Critical temperature with fixed distance between the layers versus the reciprocal of the magnetic length. Solid line 
A = 0; dashed line - |A| = 1. 



will move, which will result in energy dissipation [43j whose magnitude dissipation is proportional to the squared 
tunneling amplitude. We note that a configuration where the vortices are stationary and dissipation is absent can 
be realized fS] but such configurations cannot be used to transfer energy from source to load ^45*1. Consequently, 
systems with a negligibly small interlayer tunneling amplitude are needed. In a double quantum well the tunneling is 
quite large because of the relative narrowness and smallness of the barrier. 

The second inequality is the condition for using the so-called lowest Landau level approximation (where transitions to 
inactive levels are neglected). In this case the electrons and holes in the active level are well-determined quasiparticles. 
In systems with a quadratic dispersion law for charge carries the splitting between the Landau levels ujo — heB/m^,c 
(to* is the effective mass of the carriers) is directly proportional to the magnetic field and the Coulomb energy 
Ec ~ /el — [e^ I e) \J eB I he is proportional to the square root of the magnetic field. Consequently, the inequality 
Ec < ujQ holds only in quite strong fields (B > 10 T for GaAs). The magnetic length is small, i.e. the critical value 
of d is small. Correspondingly, large values of d cannot be used to suppress interlayer tunneling. 

For a bilayer graphene system ujq = \/2hvF/£, i.e. the condition Ec < uJc is equivalent to the condition e > 

/{y/2HvF) ~ 1.5. For an appropriate choice of e the condition is satisfied for arbitrary magnetic fields and large 
values of d can be used (the ratio d/£ can be made less than 1 by an appropriate choice of the magnetic length). Since 
graphene is a monolayer and not a sufficiently wide quantum well, the interlayer tunneling amplitude will be much 
less than for the same values of d. 

We note that the inequality Ec <^ ujq is not so fundamental for a graphene system as for a bilayer electronic system 
in a GaAs heterostructure. In the second case, in the limit of weak magnetic fields we switch to the case of two layers 
with carriers of the same t y pe , and in the first case (in a gate field) we arrive at an electron-hole graphene system 
in a zero magnetic field [2J,[25, 2^, 27 1. For e — 1 and when a large number of positive Landau levels are filled in a 
single layer (and the same number of empty negative layers in the other layer) the approach used in the present work, 
strictly speaking, is not applicable. Nonetheless, the pairing effect apparently remains. The question of the transition 
temperature in this case remains open for the time being. 

In closing, we shall discuss the question of which maximum superconducting transition temperatures can be reached 
in a bilayer quantum Hall system based on graphene. Let us assume that we have created a system with a prescribed 
value of d and we start to vary the magnetic field (at the same time adjusting the gate field so as to remain in the 
regime with zero imbalance of the active component). To determine the maximum transition temperature we shall 
represent the plot shown in Fig. [7] in a somewhat different form; specifically, we shall choose the quantity e^/ed as 
the unit of temperature. The results are displayed in Fig. [H] As one can see from the dependences presented, the 
maximum temperature is Tc ~ O.Ole^/ed and is reached for d/£ — 0.5-;- 1 in the case where carriers in the zeroth level 
participate in the pairing. For d ~ 50A and e — 3.9 the maximum temperature Tc ~ 8 K. The required magnetic 
fields S > 25 T. A transition temperature of the order of 1 K can be attained with d w 400A (for the same e) in 
sufficiently strong fields B > 0.4 T. 



APPENDIX. Spectrum of collective oscillations with arbitrary direction of the wave vector 



The problem of the collective modes in a bilayer quantum system with interlayer phase coherence is examined in 
Refs. [a, 0, [l3| and [s^. In Ref. [sl] the approach of Ref. [l3| (based on the quantization of the energy of t he p hase 
fluctuations and local imbalance) is extended to the case of an arbitrary average imbalance. It is shown in [3a | that 
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the spectrum possesses electron- hole symmetry and becomes a Bogolyubov spectrum in the low-density limit. The 
case where the direction of the wave vector coincides with the direction of the phase gradient has been examined in 
Refs. [l3] and [sl]. In Refs. [1] and [§| the functional integral approach is used to find the spectrum. In Ref. [1] the 
case of zero imbalance was examined and the spectrum was found for an arbitrary direction of the wave vector. The 
general case is studied in Ref. In the present Appendix we extend the approach used in Ref. [sl] to the case of 
an arbitrary direction of the wave vector. 

Using the analogy with a quantum ferromagnet [6] we shall examine the Fourier components of the pseudospin 
density operators 

m(q)= ^(/ci|e-^'in^2)<fc,Tya^-,fe (A.l) 

where cTy are Pauli matrices. The state with nonzero phase gradient can be interpreted as a helical state with 
rriz — cos 'm-x ~ sin^o cos(Qr), ruy = sin 0o sin(Qr). In such a state the Fourier series expansion of the averages 
{^\rhx(y){ci)\'i>) contains terms with q = ±Q. 
The wave function 



I*) = n i<k + e''^' i<k) |0> (A.2) 

describes the helical state with Q = Qx, where i is a unit vector along the x axis. To describe the helical state with 
arbitrary Q — (Qx, Qy) the wave function must be chosen in the form 



*> = n i< + e^""^' f 4fc-Q„ |0) • (A-3) 



It is easily verified that the energy is independent of the direction of Q: in the state (jA.3[) . just as in the state (|A.2[) . 
the expression for the energy has the form (|3ip . Indeed, the transition from (|A.2p to (jA.3[) is simply a rotation of the 
coordinate system, and the energy should not change under such a transformation. Wc shall examine the fluctuations 
of 9 and tp along the x axis in a new coordinate system. To take account of such fluctuations we write the wave 
function in the form 

I*) = n {^^' Y<fc + ^^'^^'^^^'^ y <^-Q„) (A-4) 

A calculation of the energy in the state (jA.4p gives 

E = Eo - eVY,cos9x + ^ J2 { i^Qy (X - X') - Fs{X - X')] cos 9x cos 9x' 

-Fd,q„(X - X')sm9x sinOx' cos(Q,(X - X') + ^x->Px')}, (A.5) 



where 



F^.qM) = Ye r /^I'x^n, e-'^v/^WF-4^,-.,.Q„^^ (^fP+xyp) . (A.6) 



2£ J-ao y/q^+Xyi'i 

Performing the expansion (|A.6[) with respect to small fluctuations ipx and mz{X) — cos9x — cos^o and switching 



to the Fourier components 



rnM-^Y.^AX)e'^^'', ^(g) = ^ ^ (A.7) 



^ Y ^ X 

we obtain the following expression for the energy of the fluctuations: 
S 



Efi = ^ ^m:,{~q)K.zz{q)m^{q) + (p{-q)lC^^{q)(p{q) ~ {imz{~q)ICz^{q)ip{q) + c.c.) , (A 

9 
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where 



/C.. (g) = Hq^ (q) - Ts (g) + (Q) + L^d (Q) 



ICz^iq) = cos 6*0 



^A\qx + Q\) + TuHqx - Q) 



cot^ 6*0, 



■^D(ki + Q|)-^D(|gi-Q|) 



^d(Q) - 



■^D(ki + Q|)+^D(|g£-Q|) 



.2!ii 1 - e-'*!?! cos{QyqP) 



(A.9) 

(A.IO) 
(A.ll) 

(A.12) 



and the functions jrs,D(<z) are determined by Eqs. (|28p and ((29)) . To quantize the energy (|A.8[) we take account of 
the fact that mz and ip are canonically conjugate and the commutator of the operators of these quantities equals 



[rhz{q),(p{q')] 



9 -2^^' A 



Expressing the operators rhz{q) and (p{q) in terms of Bose creation and annihilation operators 

„ , , ^ /27r£2 .1 /27r^2 



(A.13) 



(A.14) 



replacing in Eq. (jA.SP the quantities and (p by the operators (jA.14p . and requiring the terms containing two 
creation or two annihilation operators to vanish (which makes it possible to determine Aq), we arrive at the equation 



where 



Hn = Y,E{q) fb+{q)biq) + l 



E{q) = 2 (^^K^M^^M + ^zM^ 



(A.15) 



(A.16) 



The spectrum (|A.16|) of the collective modes depends on g, Q, and the angle between the x axis (which is chosen in 
the direction of the wave vector q) and the vector Q. Rotating the coordinate system to that the x axis is oriented 
along Q, we obtain the spectrum of the collective modes for an arbitrary direction of the wave vector: 



E{ci) - Veq,Q(eq,Q + 27q,Q sin^ Bo) + cos 6*0 (J^D(|q + Qx\) - J^odq - Qx\)) , 
where the quantities Cq^g and 7q,Q are determined by the expressions (j49[) and (jSOp . respectively. 



(A.17) 
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